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Abstract 

A model for the dynamics of actin filament ends along the leading edge of the 
lamellipodium is analyzed. It contains accounts of nucleation by branching, of de¬ 
activation by capping, and of lateral flow along the leading edge by polymerization. 

A nonlinearity arises from a Michaelis-Menten type modeling of the branching pro¬ 
cess. For branching rates large enough compared to capping rates, the existence 
and stability of nontrivial steady states is investigated. The main result is exponen¬ 
tial convergence to nontrivial steady states, proven by investigating the decay of an 
appropriate Lyapunov functional. 
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1 Introduction 

The lamellipodium is a thin protrusion, developing when biological cells spread on flat 
surfaces. It is supported by a roughly two-dimensional meshwork of protein hlaments, 
created by polymerization of actin [6] . In steadily protruding lamellipodia, the meshwork 
exhibits two dominant directions, approximately symmetric to the leading edge of the 
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Figure 1: Lateral flow. Solid and dashed lines represent the present and, respectively, a 
future state of hlaments of the leading edge (drawing courtesy of J. Vic Small). 


lamellipodium [S], and can thus be approximated by two distinct families of hlaments. 
The meshwork is a very dynamic structure, driven by the polymerization of the hlaments 
abutting the leading edge, but also by the nucleation of new hlaments via branching away 
from old hlaments, close to their growing ends. This is responsible for the two-direction 
structure with the angle between the two families approximately equal to the branching 
angle. Finally, capping of hlaments plays a role, whence hlaments become deactivated, 
stop to polymerize, and subsequentially loose contact to the leading edge. 

This work deals with a model for the dynamics of hlament ends along the leading 
edge. New hlament ends are produced by branching, they disappear from the leading 
edge by capping, and they move along the leading edge by what is called lateral flow, a 
consequence of polymerization and of the inclination of hlaments relative to the leading 
edge [S] (see Fig. [^. Under the idealizing assumption of a constant angle between 
hlaments and leading edge and of a constant polymerization speed, the speed of lateral 
how along the leading edge is constant. In reality, however, this cannot be expected, since 
the polymerization speed is subject to various inhuences such as chemical signaling and 
mechanical restrictions due to the varying geometry of the leading edge, where the latter 
will also lead to varying angles between hlaments and leading edge. 

A complete model therefore needs to describe the positions of hlaments and of the 
leading edge. The authors have been involved in the formulation of such a modeling 
framework, the Filament Based Lamellipodium Model (FBLM) [3ll3], including accounts 
for hlament bending, cross-linking, and adhesion to the substrate, as well as a number of 
other relevant mechanisms. The present work is concerned with a submodel describing 
branching and capping, and where the lateral how speed along the leading edge will be 
considered as given. As a further model simplihcation, the lateral how speed of both 
hlament families will be assumed equal at each point on the leading edge. Concerning 
the geometry, two diherent situations will be considered: For cells surrounded by a lamel¬ 
lipodium, the leading edge is described as a one-dimensional interval with a periodicity 
assumption, where the two ends are identihed. On the other hand, for cells like the mov¬ 
ing hsh keratocyte with a crescent-like shape and a lamellipodium only along the outer 
rim, the leading edge is represented by an interval with zero lateral inhow. 

The branching process requires the Arp2/3 protein complex connecting the old and 
the new hlament at the branch point [7j . The assumptions that the availability of Arp2 /3 
is limiting and that the Arp2/3 dynamics is fast compared to the branching dynamics 
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results in a Michaelis-Menten type model, similar to the one already formulated in [2]. 
Capping is described as a simple Poisson process. 

The main question of this work is: Does the mathematical model describe a stable 
distribution of hlament ends? The answer is a conditional yes with the rather obvious 
condition that the branching rate has to be big enough compared to the capping rate. 
Otherwise the hlament end population dies out. 

The rest of this paper is structured as follows: In Section the derivation of the model 
is described. This has already been explained in the context of the full FBLM in [3], but 
it is also included here for the sake of completeness. In Section]^ an existence, uniqueness, 
and boundedness result is proven. It is also shown that it is enough to initially have a 
small amount of hlament ends of only one family, to make the densities of both families 
positive everywhere within hnite time. The short Section is concerned with the proof 
of the simple result that the end distributions converge to zero, when the branching rate 
is too small compared to the capping rate. In Section existence results for non-trivial 
stationary states are proven. There are several kinds of results. First, it is shown that a 
transcritical bifurcation away from the zero steady state occurs, when the ratio between 
the branching rate and the capping rate exceeds a critical value. This local result is 
extended in two special situations: In the case of the periodic leading edge, existence of 
a nontrivial steady state is proven also far from the bifurcation point, if the lateral flow 
speed is almost constant. The same result holds for the mathematically more difficult case 
of a leading edge with zero lateral inflow, if the lateral flow speed is constant. Finally in 
Section it is shown that for every nonvanishing initial distribution the solution converges 
exponentially to the nontrivial steady state, if it exists. 


2 Derivation and Nondimensionalisation of the Model 

A model very similar to the one considered here has been formulated in [2]. We shall 
follow the derivation given in [3] in the framework of the FBLM. 

The leading edge is assumed as a (potentially closed) rectihable curve of length L, 
parametrized by arclength x G [0,L]. We distinguish between two families of hlaments, 
those pointing to the right with number density of ends u{x,t) and those pointing to 
the left with density v{x,t). By lateral flow the right-pointing hlament ends are moved 
to the right and the left-pointing hlament ends to the left, both with the prescribed 
position dependent speed c{x) >0. It is a simplifying assumption that the speed is time 
independent and the same for both families. The density of Arp2/3 is denoted by a{x, t). 
It is assumed to be recruited from the cytoplasm to the leading edge with a constant rate 
Crec with the opposite reaction working towards the equilibrium value Oq. Furthermore 
Arp2/3 is consumed by branching events, where hlament ends of one family create ends 
of the other with rate constant Kbr- Arp2/3 molecules at the leading edge are assumed 
immobile for simplicity. Finally, the rate constant for the deactivating capping reaction 
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is denoted by Reap- These assumptions lead to the system 


dtu + d^{c{x)u) 
dtv - dx{c{x)v) 


dta 


a 

Kbr ^ ^cap^ 5 

ao 

a 

Kbr U l^ca.pV , 

ao 



l^hr - {U + V) 

Oq 


We introduce the scaling 

X ^ X L, t —)■-, c{x) ^ L Kcapc{x), {u, v) ^ (u-^, v-^] , a —)■ aoo, 

^cap \ ^br ^br J 

and the dimensionless parameters 

^br ^cap^O 

a := -, e = —— , 

^cap ^rec 

where a, the ratio between the branching and the capping rates, is assumed of moderate 
size, whereas £, the ratio between the characteristic time for Arp2/3 and that of the 
capping and branching processes will be assumed as small. Whereas the hrst assumption 
is justihed and actually necessary, as our analysis will show, the smallness of £ has, to 
the knowledge of the authors, not been verihed experimentally. The nondimensionalized 
system has the form 


dtU + dx{c{x)u) = aav — u , 
dtV — dx{c{x)v) = aau — v , 

edttt = 1 —a(l + M + n). 


The last step in the model derivation is to pass to the quasistationary limit £ —)■ 0 in 
the equation for a, which is analogous to the derivation of Michaelis-Menten kinetics. 
Elimination of a from the resulting system gives 


dtU + dx {c{x)u) 
dtV - dx (c(x)v) 


a V 

1 + U + V 
au 

1 + U + V 


( 1 ) 


for X G [0,1]. Two types of boundary conditions are biologically relevant. In the case of 
a ring-shaped lamellipodium around the whole cell we assume periodic boundary condi¬ 
tions. On the other hand, if we consider only a lamellipodium at the front, it is reasonable 
to assume that no left-moving hlaments enter from the right and vice versa. These con¬ 
siderations allow to complement ([^ with one of the following sets of boundary conditions: 


(DBG) M(0,t) = 0, f(l,t) = 0, for f > 0, (2) 

(PBC) u{0,t) = v{0,t) = for f > 0. (3) 
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Throughout this paper, we will use the abbreviations {DBC) and {PBC) for Dirichlet 
Boundary Conditions and Periodic Boundary Conditions, respectively. For {PBC) we im¬ 
plicitly assume that also the lateral flow speed c{x) is periodic. To complete the dehnition 
of the problem, we pose initial conditions 

u{x^ 0) = ■Uo(2^); ^{x, 0) = ^ 0 ( 3 ^) 5 for x G [0,1] , (4) 

with given mq( 2 ^), 1 ^ 0 ( 2 ^), which are assumed to be non-negative and to satisfy the boundary 
conditions. 


3 Existence, Uniqueness, and Positivity of Solutions 


We start with a reformulation of the problem, which will be useful in most of our proofs. 
It is based on the assumption that the lateral flow speed is not only positive, but bounded 
away from zero: There exist positive constants c, c, such that 

0 < c < c < c in [0,1]. (5) 


An average inverse speed can then be dehned by 

1 da; 

C Jq c{x) ■ 

Now we introduce the transformation {x, u, v) GG {X, U, V) by 

X = C e [0,1] , cu = CU, cv = CV. (6) 

ao 

The transformed version of Q reads 

dtU + CdxU = Ruy - U , 

dtV-CdxV = Rv,u-V, (7) 


with 

R IX 

l+fi{X){U(X,t)+V{X,t))’ 

and (3 = C/c, bounded from above and below by 

0 < /? := Cjc < (3 < CIc =: (3 in [0,1]. (8) 

Note that (17, V) satishes the same boundary conditions (|^, ([^ as (u, v). The transformed 
initial conditions read 

17(t = 0) = 17o:=^, I/(t = 0) = Vo:=^, in [0,1] . (9) 

The reformulation has two effects: First, the solution of the equations by the method 
of characteristics is simplihed, since the transformed lateral flow speed is constant, and, 
second, linearization around the zero solution gives a problem with constant coefficients. 
Applying the method of characteristics leads to a mild formulation. 
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Definition 1. (PBC) Let {U,V) G C(T^ x [0, oo))^ (with the torus represented by the 
interval [0, 1]) satisfy 

U{X,t) = Uo{X-Ct)e-^+ [ e^-^Ru,v{X+ C{s-t),s)ds, 

Jo 

V{X,t) = Vo{X + Ct)e-^+ [ e^-^Rv,uiX-C{s-t),s)ds, (10) 

Jo 

for {X,t) G X [0, oo). Then 

{u{x,t),v{x,t)) = -^(U(X(x),t),V(X(x),t)) (11) 

c(x) 

is called a global mild solution of the problem (0, <Q, fl. 

(DBC) Let {U, V) G C([0,1] x [0, cx)))^ satisfy 

U{X,t) = Uo{X-Ct)H{X-Ct)e-^ + [ e^-^Ru,v{X + C{s - t),s) ds, (12) 

J{t-X/C) + 

V{X,t) = Vo{X + Ct)H{l - X - Ct)e-^ + [ e^-^Rv,uiX - C{s - t), s) ds, 


for {X,t) G [0,1] X [0, cxd), where H denotes the Heavy side function, 
by (11) is called a global mild solution of the problem 0, 0. 


Then {u, v) defined 


Proposition 1. Let Uq, Vq, c G C([0, 1]) satisfy ( [5|, Up, Vq > 0, and the boundary conditions 
(DBC) or (PBC). Then the problem 0, 0(~0, or, respectively, 0, 0, 0 , has a 
unigue, global mild solution {u,v), satisfying 


0 < u(x, t) <- max < maxwo, Oi \ , 0 < vix, t) <- max < max vq, a \ . (13) 

C ( [0,1] J C ( [0,1] J 

Remark 1. It is straightforward to show by differentiation of the eguations that for smooth 
initial data uo,vo G C°°(T^) for (PBC), and uq G C“((0, 1]), vq G C“([0,1)) for (DBC), 
and for smooth lateral flow speed c G C“([0,1]), the solution satisfies u,v G C“([0,1] x 
[0,oo)). Some results of the following sections are based on computations with the strong 
forms 0 or 0 o/ the differential eguations. These can be justified by uniform smooth 
approximations ofuo,vo,c, and subseguent removal of the smoothing. 


Proof. The (obvious) non-negativity and Lipschitz continuity of Ruy in terms of (f/, V) G 
[0, oo)^, as well as the bound Ru,v ^ Oi/fl will be sufficient for carrying out the proof. 

With the Lipschitz continuity, it is straightforward to show that the right hand sides 
of (10) and (12) preserve non-negativity and are contractions on C(Tx [0,T]) and, respec¬ 
tively, C([0,1] X [0,T]), for T small enough, which proves local existence. The estimate 


0 < u{x, t) < 


C 

c 



C \ 

< — max < 

sup Uo, y . 

1 c ] 

> < - max < 

maxun,« [ 

c 1 

[ [0,1] ^ J 

1 c 1 

1 [0,1] j 


and the analogous version for v allow to continue the local solution indefinitely and also 
prove ([10 . □ 
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One expects that if the initial conditions '^^ 0 ( 2 ^)) are positive (except at the 

boundaries for (DBC)), the same holds for the solution {u{x,t),v{x,t)) for all f > 0. In 
fact a much stronger result is true: It is enough to have positivity of the initial data on 
some interval for only one family; after hnite time both families will be positive everywhere 
(except at the boundaries for (DBC)). The reason for this is that, although the initial 
mass is reduced by capping while it is transported across the domain, it remains positive. 
This mass, however, will trigger the creation of new hlaments of the other family through 
the branching term. This new mass will be transported in the opposite direction and will 
itself cause the creation of mass of the hrst family. As long as the lateral flow speed c{x) 
is bounded from below, this process happens in hnite time. 

Proposition 2. Let the assumptions of Proposition^ hold and let the initial data satisfy 

{uq + vq) dx > 0 . (14) 

Then, for every T > 2 dx/c(x), 

for (PBC): u,v > d in [0,1] x [T, 00 ), 

for (DBC): m > 0 in (0,1] x [T, cx)), n > 0 m [0,1) x [T, cx)). 




Figure 2: Illustration of the proof of Proposition!^ 


Proof. The equivalent result for {U, V) will be proved. Assumption (14) and the continuity 
of the initial data imply the existence of an interval [Xi, X 2 ] C [0, 1] of positive length and 
of m > 0, such that one of the initial densities, w.l.o.g. Fq, satishes U^Pniin [Xi,X 2 ]. 

The mild formulations given in Def. [^show that, whenever a 17-characteristic (with 
velocity C) passes through a region where V is positive, U is positive along this charac¬ 
teristic in and after (timewise) this region, and vice versa {U ■H- V, C ^ ~C). Also U is 
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positive everywhere along a 17-characteristic after a point on the characteristic, where it 
is positive (again the same for V). 

These observations reduce the proof to a geometric problem (see Fig. [^, where we 
alternate between using the equations for U and for V. The hrst step is the observation 
that by the above property of the initial data, U is positive in the strip iSq, defined by 
Xi < X — Ct < X 2 (light gray shading in Fig. [^. This implies, as the second step, 
that V is positive there and along all 17-characteristics starting in this strip, i.e. also 
in the region 5i bounded by X = 0, X — Ct = Xi, X + Ct = Xi, X -|- Cf = 2 — Xi 
(gray shading in Fig. |^. The third step is to draw again U characteristics, now starting 
in iSi, which adds a triangle S 2 (dark gray shading in Fig. above iSi, where U > 0. 
Straightforward continuation shows that for t > T* = (1 — Xi)/C -|- 1/(7, U and 17 are 
positive for 0 < X <1. Positivity on the boundary, except for 17(0, f) = 17(1, f) = 0 in the 
(DBG) case, is achieved at any time after T*, implying that the result of the proposition 
holds with any T>2/(7>T*. □ 


4 When Capping Exceeds Branching 


The dimensionless parameter a is the ratio between the branching and the capping rate. 
If it is too small, it can be expected that u and v tend to zero as t —)■ 00 . Note that, 
on the other hand. Proposition holds for all a > 0. This means that even if we start 
with initial conditions, which are zero everywhere and are only positive for one family in 
a small interval, the solutions hrst become positive everywhere before they decay to zero. 

Proposition 3. Let the assumptions of Proposition and 0 < a < 1 hold. The the 
solution (u, v) of (§, (§, (§ or Q, (§, (g satisfy 

ll'“(■)^)lln2((o,l)) + ll'^(■)^)llnh(o,l)) - (lko||L2((o,i)) + lko||i,2((o,i))) • 


Proof. We observe that 




2auv 
1 + U + V 



dx 


0 

—c{l)u{l,t) — c(0)n(0,t) 


(PBC), 
(DBG). 


The non-positivity of the boundary term and the estimate 

2auv o o {a — l){u‘^ + — a{u — vy — {u + v){u‘^ + 0 “^) 


1 + U + V 


2 2 
— U — V = 


< (a — 1) {u^ + 


1 + U + V 


hnish the proof. □ 

The results of the following sections show that the bound on a is sharp for (PBG), 
but not for (DBG), where decay to zero can be expected also for 1 < a < oq with the 
bifurcation value oq- 





5 Existence of Nontrivial Steady States 


Bifurcation from the zero solution 

A bifurcation value ao > 1 for the ratio a between the branching rate and the capping 
rate will be computed, where a nontrivial steady state branches off the zero solution. The 
above inequality is derived easily by repeating the computation in the proof of Proposition 
[^for the steady state problem, proving that for a < 1 it has only the zero solution. 

We shall work again with the variables {X, U, V), where the linearization around the 
zero solution has constant coefficients. We rewrite ([^ as 



(15) 


with the linearized operator 



aV-U- CdxU 
aU -V + CdxV 


The bifurcation point will be the smallest value ckq > 1, where £^ 0 ) subject to the bound¬ 


ary conditions (|^ or (|^, has a nontrivial null space. 

By inspection it is obvious that for (PBC) ao = 1 holds with ([/, V) := (1,1) G Af{Ci). 
Straightforward computations show that the null space is one-dimensional. 

For (DBG) another straightforward, although a little longer computation gives that 
bifurcation values are of the form a = l/ \ cos6|, where b solves 


(16) 


Cb -|- tan 6 = 0 . 


We denote by bo G (7r/2,7r) the smallest positive solution, set cxq = l/|cos6o|, and note 
that 



Note that [/ > 0 in (0,1] and fo > 0 in [0,1). Again the null space is one-dimensional. 

For both types of boundary conditions, inhnite increasing sequences of bifurcation 
values exist. However, for all bifurcation values larger than ckq, the null spaces, and 
therefore the bifurcating solutions, consist of functions with changing signs, which are 
irrelevant for our application and have no chance to be the long-time limit of non-negative 
solutions. 

The normal form reduction close to the bifurcation point is derived by choosing values 
of a close to ao and by making the ansatz 



(18) 


For the linearized operator and its formal adjoint with respect to the scalar product in 
L^((0,1))^, the symmetry property 
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holds. Thus, the scalar product of (15), written in the form 

a^{U + V) 


a 


with (K, U) gives 


= £„„(17,l/) + 


a 


tto 


1 + I3{U + V) 


— / {UV + VU) dX = 

d.t ./n 


Ol — Qq — 


al3{U + V) 

1 + P{U + V) 


{VV + uil) dX . 


Substitution of the ansatz (18) leads to 

dB 


dt 


= (a — ao)B{Ki — 0^2) + 0 {{a — ao) ), 


(19) 


with 


fhu^ + V^)dX 

= >0, 


aoj;m + V){U^ + V^)dX ^ ^ 

K-2 — - „i ~ ~- > U . 


2f‘UVdX 2f„‘UVd.X 

This is the normal form of the transcritical bifurcation. It indicates that for a < ao the 
trivial steady state is stable, whereas for a > ao, stability is transferred to the bifurcating 
steady state 


1 / 


= {a- ao) — 
^2 


1 / 


+ 0{{a — ao) ) 


( 20 ) 


This result is of course only formal. For the steady state bifurcation, however, our compu¬ 
tations, in particular ki,K 2 7 ^ 0, verify the conditions of the Crandall-Rabinowitz theory 
[1] for bifurcations from a simple eigenvalue. Without going into further detail, we shall 
state the result below. We do not attempt to rigorously justify the dynamic normal form 
reduction (19), since the stability of the bifurcating state will be proved in the following 
section, anyway. 

Proposition 4. Let ^^hold and let ao be defined as above. Then there exists ai > ao, 
such that the system^^ with ([^ or ([^ has a smooth (with respect to a) branch of 
nontrivial solutions {u,v) : [ao,ai) —?• C([0,1])^ with {u,v) 1 _ = (0,0) and with 


u 

da Vh 


[X 




(U, V) c 


I Q;=Q;o 

dy 


'0 


c{y) 


a=ao K,2C 

In a neighborhood of the point {u,v,a) = (0, 0, ao) in C([0,1])^ x M no other stationary 
solutions besides the trivial solution and the solutions on the nontrivial branch exist. 

Lemma 1. Let the assumptions of Proposition^hold. Then for a —ao > 0 small enough, 

for (BBC): u,v > 0 in [0,1] , 

for (DBC): h > 0 in (0,1], h > 0 in [0,1) . 

Proof. Proposition justihes (20), with 0{{a — ao)^) to be understood uniformly in 
X G [0,1]. This immediately implies the result, except for the case of (DBC), where the 
behavior of U close to X = 0 and of V close to X = 1 has to be examined. However, 

dC 


r(o) = o, -^(')) = 0r(o)>o, C(i) = o, 


a . 


d.Y(l) = -X(l)< 0 ' 


implies positivity of U near X = 0 and of V near X = 1, completing the proof. 


□ 
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Periodic boundary conditions — almost constant lateral flow speed 

In the case of periodic boundary conditions and of a constant lateral flow speed, the 
bifurcating branch of nontrivial solutions is given by 


u{x)\ _ a — 1 /l\ 

vlx)J 2 


a > 1 . 


( 21 ) 


So it is explicit, homogeneous, and global. We make use of these properties to construct 
a global branch for almost constant lateral flow speeds, satisfying 


c(x) = Co + eci(x ), 


with a constant Cq > 0, with a smooth function Ci(x), and with a small parameter e. 
Substitution of the transformation 





1 ) 




in Q gives 


AM 
dx \Vi J 



h{x) + er{ui,Vi,x] a, e ), 


( 22 ) 


with 


M(a) 


1 A - 3a 
2 coa 1 


a + 1 \ 
3a-y ’ 


h{x) 


c[{x) 

2co Uy* ’ 


(23) 


and where r is smooth in all its variables and uniformly bounded in terms of a —)■ cxd and 
e —)■ 0 with commuting limits. It also satishes 


r(a 


1) = r(tii,!i,,i;ir) 



(24) 


These observations are the result of a lengthy but straightforward computation. The next 
step is to rewrite (22) as a hxed point problem by inverting the linear operator on the left 
hand side on the space of periodic functions. An explicit computation gives 


/ Ui{x; a, e) 
Ai(x;a,e:) 



e ^^{h + er){y) dy + 


e ^'^{h + er){y) dy 




(25) 

requiring the existence of the inverse of e ^ — I, equivalent to the invertibility of M, which 
is true for a > 1. Since M(a) also converges to an invertible matrix as a —)■ oo, the linear 
operator applied to h + sr on the right hand side is uniformly bounded in a G [a*, oo) 
with a* > 1. Obviously, for £ small enough, solutions can be constructed by contraction 
on the space ^(Tr^ x [a*, oo))^ of bounded continuous functions of the periodic variable 
X and of the parameter a. We skip the details of the proof. 
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Proposition 5. Let a* > 1 and let 

c{x) = Co + eci{x ), with cq > 0 , Ci G C^(T^). 


Then there exist Eq^C > 0 such that for every e < Eq and for every a > a*, the problem 
0 . @ with (PBC) has a smooth stationary solution {u,v), satisfying 

lltt — (a — l)(l/2 + EUi{- ; a, 0 ))||l°°([o,i]) + 11 ^ — (o — 1)(1/2 + evi{- ; a, 0))||i;,oo([o,i]) < e'^aC , 


where (mi(x; a, 0), a, 0)) E Cb(T^ x [a*,cx)))^ is given by (25) withE = 0. 

Remark 2. Note the uniformity of the result in terms of a E [a*,oo). In particular 
the upper bound Cq for e and the constant C in the error estimate are independent from 
a -E oo. There is some subtility to the situation as a —)■ 1, since the matrix —1)“^ 

blows up in this limit. However the limiting right hand sides being proportional to the 
vector (1,1) (see (23), (24) ) satisfy the solvability conditions for a = 1, so that a bounded 
passage to the limit can be expected. We do not carry out this limit in detail, since it is 
roughly eguivalent to the bifurcation analysis above. 


We conclude with an example: Setting 


yields 


Mi(x; a, 0) = 

Vi{x; a, 0) = 


c{x) = 1 ^ cos(27ra:) 


Stt (2Q;7r cos(27ra;)) — (a — 1) sin(27ra;) 
8 (a — 1 + 27r‘^a) 

Stt (2Q;7r cos(27ra;)) + (a — 1) sin(27rx) 
8 (a — 1 + 2Ti‘^a) 


(26) 


(27) 


Figure shows the approximate steady state solution (a — l)(l/2 + eui{x; a,0), 1/2 + 
EVi{x] a,0)) together with a numerical solution computed as a steady state of the time 
dependent problem. The lateral flow in the example corresponds qualitatively to a simula¬ 
tion with the full Filament Based Lamellipodium Model [3], where the cell is moving and 
the front is located at x = 0 (identihed with x = 1) and the back at x = 0.5. The lateral 
flow speed c(x) has its minimum at the back and its maximum at the front. The steady 
state distributions of both hlament families have their maximum at the back, due to the 
accumulation of hlaments by the lateral flow. The maximum of the right-moving family 
is shifted slightly to the left as compared to the cell rear and vice versa for the left-moving 
filaments. The reason for this seemingly counter-intuitive result is that filaments of the 
right-moving family are produced by left-moving filaments, which shifts the maximum of 
the right-moving family to the left. 


Dirichlet boundary conditions — constant lateral flow speed 

In terms of the new unknowns p = u + v,q = u — v, the stationary version of ([^ with 
constant lateral flow speed c can be written as 


cp' = —q 



cq = p 


a 

l+p 



(28) 
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Figure 3: Approximate and numerical solution for a = 10 and £ = 0.1 with c{x) given 
by (26). The end density of the right moving hlaments u is depicted in blue (dashed) 
and that of the left moving hlaments v (solid) in red. Thin lines are the asymptotic 
approximations, thick lines are the numerical solution of the time dependent problem 
after t = 20 and the left |/-axis applies. In green (dotted) the lateral how c{x) is depicted, 
the values are in relation to the right y-axis. 



The boundary conditions (DBG) translate to 

p(0) + g(0) = 0, p{l)-q{l) = 0. 

The bifurcating solutions constructed above have the symmetry 

u{x) = h(l — x) i — > p{x) = p{l — x ), q{x) = —g(l — x ), 

and we shall look for solutions with this property. This allows to reduce the problem to 
the interval [0,1/2] with the boundary conditions 

p(0) + g(0)=0, q{l/2) = 0. 


Viewing (28) as a dynamical system and assuming a > 1, it has the critical points 
(p, q) = (0, 0), which is a center, and the saddle (p, q) = [a — 1, 0). We look for solutions 
following trajectories in the region bounded by p + g = 0, by g = 0, and by the stable 
manifold of the saddle. The system has the hrst integral 


Eq = q^ + E(p, a) , with E(p, a) = —p^ + 4ap — 4a(a + 1) In ( 1 H- - — 

V a + 1 

which has been used for drawing the trajectories in Fig. 

The above mentioned region corresponds to 

Eq G [0, E*{a)] , with E*{a) = E{a — l,a) = (3a + l)(a — 1) — 4Q;(a +1) In 


( 29 ) 


2q; 

cr T 1 
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Figure 4: This figure shows the trajectories in {p, gj-phase space of the stationary equation 
with constant c and {DBC). The orange, diagonal lines represent the boundary conditions 
q = —p at a: = 0 and q = p at x = 1. Dashed curves represent (29) for different values of 
Eq and the blue lines the corresponding solutions curves fulfilling p(0) + g(0) = 0. The 
red line is desired the solution satisfying both boundary conditions. 



where = 0 corresponds to the origin and Eo = E*{a) to the stable manifold of the 
saddle. Among the trajectories connecting the segments p + q = 0 and g = 0, we need 
one which takes ’time’ 1/2. To find an appropriate value of Eq, we compute g < 0 from 
(29), substitute it in the first equation in (28), and integrate. This produces an equation 
for Eq: 

fPiiEo,a) ^ I 

I{Eo,a) = c - ^ , (30) 

Jpo{Eo,a) (1 + a/{l+p))^^Eo - E{p,a) 2 

where p(0) = po{Eo, a) and p(l/2) = pi{Eo, a) are the unique solutions of 

Eo = E{po,a)+ pI, Eo = E{pi,a). 

Information about the range of I{-,a) is needed. For studying the limit as Eq —)■ 0, the 
Taylor expansion E{p, a) = + 0{p^) and its consequence po ~ ^yEo{a + 1)/(2q;), 

Pi ~ ^yEo{a + l)/{a-l), imply 


J(0, a) = 


-\/a^ — 1 


71 

-arcsm 

2 


a 


2a 


As a consistency check, it is easily verified that the unique solution of the equation 
/(0,Q!) = 1/2 is the bifurcation value ao defined by (16), (17). The uniqueness follows 
since /(O, a) is strictly decreasing as a function of a, i.e. 

/(O, a) < 1/2 iff a > ao . 


On the other hand, by pi{E*{a), a) = a — 1 and dpE{a — l, a) the integrand in (30) looses 
its integrability as Eq —)■ E*{a), implying 

lim I{Eq, a) = oo , 

Eq^E* (a) 
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which is actually obvious, since along the stable manifold of the saddle the critical point 
and (therefore the line (? = 0) cannot be reached in finite ’time’. These results and the 
continuity of I{Eq, a) complete the proof of the following proposition. 

Proposition 6. Let a > ao as defined by (0, 0 and let c{x) = const. Then the 
equation ([^ together with (DEC) has a non-trivial stationary solution. 

Remark 3. Uniqueness of the steady state is equivalent to strict monotonicity of I{Eo, a) 
as a function of Eq. This will be proved indirectly by the decay result of the next section. 


6 Stability of Nontrivial Steady States 


The main result of this section is that whenever a nontrivial steady state, as constructed 
(for certain cases) in the previous section, exists then it attracts all nontrivial solutions. 
We therefore make the assumption that there exists a stationary solution {u{x),v{x)) of 
([^, satisfying 

either (PBC), (u, u) G C^(T^)^, and m < u{x),v{x) < M, (31) 

or (DBG), (h, h) G C^([0,1])^ , and m \ ^ < M, (32) 

X 1 — X 


for X G [0,1], with positive constants m,M. 

For proving exponential convergence to nontrivial steady states, we need to strengthen 
the result of Proposition Note that the following result is sharp in terms of the values 
of the parameter a, proving uniform-in-time bounds away from zero for all parameter 
values above the bifurcation value, where the zero solution looses its stability. 


Lemma 2. Let the assumptions of Proposition hold and let a > a^, where ao is the 
bifurcation value as defined in the preceding section. Then for (PBC) there exist T,m > 0 
such that 

u{x, t),v{x, t) > m for {x, t) G [0,1] x [T, oo). 

For (DEC) there exist T, a > 0 such that 


u{x, t) > a sin 



v{x, t) > a sin 



, for {x, t) G [0,1] X [T, oo) . 


Furthermore, again for (DEC), there exists a constant M > 0 such that 


u{x, t) < Mx , v{x, t) < M(1 — x ), for (x, t) G [0,1] x [T, cx)). 
Proof. (PBC): As a consequence of Proposition]^ there exist mo,T > 0, such that 

u{x, T),v{x, T) > mo > 0 for x G [0,1] . 


We choose 


M 


c , f a — 11 
-mm|mo,—j 


> 0 , 


15 







implying 


U{X,T),V{X,T) > M for Xe [0,1], 

Assume the bound holds for V for later times and consider the equation along character¬ 
istics for U: 


U = 


aV 


U > 


aM 


U. 


1+/3{U + V) -1 + I3{U + M) 

The right hand side is nonnegative for U = M with the consequence U > M. Analogously 
for V. This completes the proof with m = MC/c. 

(DBC): By the result of Proposition there exist M > 0, Xq E (0,1), t 2 > ti > 0, such 
that V > M in [0,Xo] x [^ 1 ,^ 2 ]- Along [/-characteristics staying inside this rectangle, 

aM aM 


U > 


U > 


2{l + 2/3M) ’ 


l + /?([/ + M) 

holds, as long as, furthermore, 

TT r 71 ^ ctAf , 

U < mm < M, —-=—- )■ . 

\ ’2(l + 2/?M)j 

For characteristics starting at X = 0, t = t* E [ti,t 2 ), this implies 

aM . aM X 

-(*-*•) = 


until 


U > 


U = min < M, 


2(1 + 2(3M) 
aM 


2{l + 2ftM) C ’ 
min{l,t 2 - t*,Xo/C} 


2(1 + 2/3M) 

As a conclusion, there is a time T E (ti, [ 2 ] such that U (X, T) (and therefore also u{x, T)) 
increases (as a function of X) at least linearly away from X = 0 up to a certain point, 
after which, by Proposition it is bounded from below by a positive constant. This 
implies the existence of Aq > 0, such that 

U{X, T) > AoU{X ), V^(X, T) > AoV{X ), 

with [/, V dehned by (16), ( [l7| ), and the second inequality is proven analogously. Now we 
choose 

a — oq 


A := min < Aq, 


2aol3 


and assume ld(X, t) > AV{X) for (X, f) E [0,1] x [T, cx)), in the equation for U. This 
implies 


c2g> 


aAV 


-u, 


t>T. 


dX - 1 + (3{U + AV) 

We claim that U{X) = AU{X) is a subsolution for t > T, which follows from AU{X) < 
f/(X, T) and 

d{AU) 


C- 


dX 


= CA 60 cos( 6 oX) = — A tan 60 cos( 6 oX) 

= ao[A(sin bo cos{boX) — sin( 6 oX) cos [q)] — A sin( 6 oX) 


aAV 


1 -h (tt/tto — 1 ) 


- AU < 


aAV^ 

1 + 2,5A 


- AU < 


aAV 


1 + I3{AU + AV) 


-AU. 
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On the other hand, it can be proved analogously that U > AU implies V > AV. This 
completes the proof of the lower bound with a := AC/c. 

For the proof of the upper bound we observe that along characteristics 

dx - 

implying at most linear growth of 1/ in terms of X, and therefore also of u in terms of x, 
with the analogous argument for v. □ 


The convergence analysis is based on the Lyapunov functional 


'H[u,v] '■= ^ [ c + ^{v — vY') dx . 

^ ./n V / 


(33) 


Note that under the assumptions of Lemmaand with (31) or (32), the Lyapunov func¬ 
tional is well dehned along the solution of (^ for t > T, since for (DBG) the integrand 
can be continuously extended by the value zero to a; = 0,1. For the computation of the 
time derivative of 'H[M,n] along solutions, the computation 


c- {u — u) dx{cu) dx 


u 


= - -a 

./n CU 


^cv {cu — cuY 


dx— c—{u — u) dx{cu) dx 


u 


{cu — cuY udx{cv) — vdx{cu) 


cu^ 


dx— c—{u — u) dx{cu) dx 


u 


{u — uY 
2v? 


= / C 


{udx{cv) — vdx{cu)) dx — c — {u — u) dx{cu) dx 


u 


will be used, where the derivatives on the right hand side can be eliminated by using the 
stationary equations. The integration by parts leading to the third line does not produce 
any boundary terms since v{u — uY/u is periodic for (PBC) and the continuous extension 
vanishes at x = 0,1 for (DBG). The analogous computation for the second part of ?{[«, x] 
leads to 


_d . 1 ^ / cJ(n,x,x,x) 

df 2 Jq {1 + u + v){1 + u + v) 


dx , 


with 


J{u,v,u,v) : = 


-\ 2 
u — u 


u 


-\ 2 
V — V 


V 


[(1 -I- M -f v) {vY -f x^) -I- 2x^x] 
[(1 -I- X + x) {vY -I- x^) -I- 2x^x] 


(34) 


U 'll U %) /_9 _9 _ 9 _ _ 9 \ 

—2—z-z— [u -f x^ -I- u^v -I- uvA . 

u V ^ ' 

Nonnegativity of J is not obvious, but actually even a coercivity property can be shown: 
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Theorem 1. Let the assumptions of Proposition^ hold, let a > ao with the bifurcation 
value ao, let a stationary solution {u,v) of ([^ exist, which satisfies either (31) or (32), 
and let 'H[u,v] be defined by (33). Then there exists a constant 7 > 0 such that 

n[u,v] (t) < v] (T), fort>T, 


with T from Lemma\^ 


Proof. The representation (34) of J suggests the notation 


u — u 


V — V 


u := 


V : = 


u 


with {u,v) e [— 1 , 00 )^. This domain can be split into three subdomains, dehned by 
■uh < 0, 0 < u/i) < 1, and, respectively, 0 < v/u < 1. 

For hh < 0, J is obviously nonnegative and satishes 


J > 
> 


(-2 I -2\ /-2 I '~2\ rv _ 2 I / -\2\ 

(u + n (m + n = —— -iu — u) + ^('y — 'yj 
^ ^ ^ ' uv \u V / 

2 (^{u — uY + ^{v — . 

\u V J 


Now we consider the second case 0 < M/h < 1. The result for the third case then follows 
by symmetry. With A = u/v, we introduce another change of variables {u,v) —)■ (A,!)) G 
(0,1] X [- 1 , 00 ), giving 


^ = A^ [(1 + M + n + {uX + v)v) {u^ + + 2mw^] 

+ (1 + u-\-v + (u\ + v)v) (ti^ + + 2h^h 

—2A [if + f + ifv + uf) 

= [if + f) (mA + v) (A^ + 1 ) (-0 + 1) 

+ (1 — A) [(1 — A) [if + + if{X^ + 1) + (1 — Xfuf + 2fv\ 

> [if + f) [v{v + !) + (!- A)^) . 


With one more set of variables, u = ufuv = ■\/v/u{u — u), v = 
the above inequality is equivalent to 


vvuv 


J > 




uv 


[vif + {u — h)^) > 2{u — vY + 


KV 


with K = inf 

[ 0 , 1 ] 


= ^/uJv{y - v), 

[if + f) V 


uv 


The results of Lemma and the assumption (31) or (32) imply that if + f and v/v 
are bounded from below and u is bounded from above, with the consequence k > 0, and 
therefore 

J > - - [u^ + v^) = - - (-{u - uY + -[v -vY) ■ 

With the common upper bound M for u,v,u, and v, the proof is completed with 

2aK 

7 = 


(1 + 2M)2(4 + k) 


□ 
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